Directed transport driven by the transverse wall vibration 
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Directed transport of overdamped Brownian particles in an asymmetrically periodic tube is in- 
vestigated in the presence of the tube wall vibration. From the Brownian dynamics simulations we 
can find that the perpendicular wall vibration can induce a net current in the longitudinal direction 
when the tube is asymmetric. The direction of the current at low frequency is opposite to that at 
high frequency. One can change the direction of the current by suitably tailoring the frequency of 
the wall vibration. 
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I. INTRODUCTION 

Stochastic transport on periodic asymmetric sub- 
strates far from equilibrium has raised widespread inter- 
est in the recent literature^. In these spatially periodic 
structures, directed motion of particles can induced by 
zero-mean nonequilibrium fluctuations and noise. These 
come from the desire to understand molecular motors^, 
nanoscale friction'^, surface smoothing^, coupled Joseph- 
son junctions^, optical ratchets and directed motion of 
laser-cooled atoms*', and mass separation and trapping 
schemes at the microscale^. The fundamental condition 
for directed transport to occur is that either the spatial 
reflection symmetry of the system is broken or fluctua- 
tions are statistically asymmetric. Several models have 
been proposed to explain this transport mechanism under 
various nonequlibrium situations. Typical examples are 
rocking ratchets-, flashing ratchets^, diffusion ratchetsi^, 
correlation ratchetaii. 

Most studies have revolved around the energy barrier. 
The nature of the barrier depends on which thermody- 
namic potential varies when passing from one well to 
the other, and its presence plays an important role in 
the dynamics of the system. However, in many trans- 
port phenomena^-, such as those taking place in micro- 
and nano-pores, zeolites, biological cells, ion channels, 
nanoporous materials and microfluidic devices etched 
with grooves and chambers, Brownian particles, instead 
of diffusing freely in the host liquid phase, undergo a 
constrained motion. In these cases the entropic barriers 
may appear when coarsening the description of a complex 
system in order to simplify its dynamics. Reguera and 
coworkersiS used the mesoscopic nonequilibrium thermo- 
dynamics theory to derive the general kinetic equation of 
the motor system and analyzed in detail the case of dif- 
fusion in a domain of irregular geometry in which the 
presence of the boundaries induces an entropy barrier 
when approaching the dynamics by a coarsening of the 
description. In their recent work?^. they studied the cur- 
rent and the diffusion of a Brownian particle moving in 
a symmetric channel with a biased external force. In our 
previous worbi^, we studied the transport of overdamped 



particles in a three dimensional tube in the presence of 
unbiased external forces and the motion of particles can 
indeed be rectified, with a sign that depends on the de- 
tails of the wall profile. In the present work, we study 
the directed transport of overdamped Brownian particles 
moving in asymmetric tube in the presence of the tube 
wall vibration. We emphasized on finding whether the 
transverse wall vibration can induce a longitudinal net 
current and how the parameters of the system affect the 
transport. 



II. MODEL AND METHODS 




FIG. 1: Schematic diagram of a tube with periodicity L. The 
vibrating wall is described by its time dependent half width 
r{x, t) = a[sin(2p) + ^ sin(^)] + b + bo sm{ujt). A is the 
asymmetric parameter of the tube shape and A = —1.0 in 
the diagram. The sohd line describes the original shape of 
the tube and the dotted line shows its vibrating shape. 



In this study, we study the transport of the particles in 
an asymmetrically periodic tube driven by the tube wall 
vibration [Fig. 1]. Since most of the molecular transport 
occurs in the overdamped regime, we can safely neglect 
inertial effects. The overdamped dynamics can be de- 
scribed by the following Langevin equations written in a 
dimensionless forroi^ii^, 
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where x, y, are the two-dimensional (2D) coordinates, 
rj is the friction coefficient of the particles, fc^ is the 
Boltzmann constant, T is the absolute temperature and 
S,x.y{t) is the Gaussian white noise with zero mean and 
correlation function: < £,i{t)^j{t ) >— 2di,jS{t — t) for 
i,j = x,y. < ... > denotes an ensemble average over the 
distribution of noise. S{t) is the Dirac delta function. 
Imposing reflecting boundary conditions in the trans- 
verse direction ensures the confinement of the dynamics 
within the tube, while periodic boundary conditions are 
enforced along the longitudinal direction for the reasons 
noted above. The vibrating wall is described by its time 
dependent half width 



,2Trx. A . Anx.^ 
r{x) = a[sm( — ) + — sm( — )] 



b + bosm{ojt), (3) 



where a is the parameter that controls the slope of the 
tube and A is the asymmetric parameter of the tube 
shape, b is the parameter that determines the half width 
at the bottleneck, bo and uj are the amplitude and fre- 
quency of the vibration, respectively. 

In this case, no general valid analytical expressions are 
possible. However, the study of these transport phenom- 
ena is in many respects equivalent to an investigation 
of geometrically constrained Brownian dynamics^^. The 
behavior of the quantities of interest have been corrob- 
orated by Brownian dynamic simulations performed by 
integration of the Langevin equation using the standard 
stochastic Euler algorithm. The average particle veloc- 
ity in the x-direction has been derived from an ensemble 
average of about 3 x 10'' trajectories according to the 
following expression. From Eq. (1-2) and the standard 
stochastic Euler algorithm, the single integration steps 
readi^: 



Vr, 



Xold 



Void 



V2fcsTAti?i, 
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)|, the modification positions oi x,y- 

coordinate are 



Xold 



^newi Void Unewi 



(6) 



if 



Une 



> \r(xnew)\, the particles will meet the wall. In 



this case, we use a simple way to simulate the movement 
of particles near the wall: When the particles meet the 
wall, it will return to its previous position and the time 
step is incremented, namely. 



Xold = Xold, Void — Void, 



(7) 



where Xoid, Void are the original positions and x^ew, Vnew 
are the new positions. R2 are two Gaussian dis- 

tributed random numbers with unit variance. At is the 
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FIG. 2: Dependence of the mean velocity v on time step At. 
(a)for different values of temperature at A = —1.0, uj — 11.0, 
a = 0.5, b = 0.9, feo = 0.3, and L = 1.0. The full symbols 
depict the present algorithm and the empty symbols indicate 
the strict reflection algorithm, (b) for different parameters at 
A = -1.0, a = 0.5, 60 = 0.3, and L = 1.0. 



integration step time. The average particle velocity along 
the x-direction. 



V = (x) = lim 
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In order to check the validity of the present Brwo- 
nian dynamics algorithm, it is necessary to compare the 
present results with a real reflecting boundary condition 
on the wall. From Fig. 2(a), we can find that the re- 
sults from the present algorithm well agree with that 
from the strict reflection algorithm for very small time 
step. Moreover, the present algorithm takes much less 
CPU time than the strict reflection one. We have also 
checked the convergence of the algorithm in a large range 
of the system parameters. From Fig. 2 we can see that 
the algorithm is convergent and the numerical results will 
not depend on the time step for very small time step. 
Therefore, our algorithm will give well approximate re- 
sults with respect to the strict reflection algorithm. In 
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order to provide the requested accuracy of the system 
dynamics time step was chosen to be smaller than 10~^. 

We must point that the extension of our scheme 
to three-dimension with rotational symmetry along the 
transport axis is possible as well. This will consume more 
computation time, but the over results remain qualita- 
tively robust. 



III. NUMERICAL RESULTS AND DISCUSSION 

Our emphasis is on finding the asymptotic mean ve- 
locity which is defined as the average of the velocity over 
the time and thermal fluctuations. In order to obtain the 
mean velocity, we carried out extensive numerical simu- 
lations based on Eq. (4-8). For simplicity we set = 1 
and fcs = 1 throughout the work. 
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FIG. 3: The mean velocity « as a function of vibrating 
frequency uj for different values of asymmetric parameter 
A = -1.0, 0.0, and 1.0 at T = 0.05, a = 0.5, b = 0.9, 6o = 0.3, 
and L = 1.0. Wi is the frequency for the left peak, loc is the 
cross-over frequency and ojr denotes the frequency for the 
right valley. 



Figure 3 shows the mean velocity v as a function of 
the vibrating frequency ui. From the Fig. 3, we can see 
that for low frequencies the current is positive for A < 0, 
zero at A = 0, and negative for A > 0. The phenomenon 
will be opposite to that for high frequencies. It is obvi- 
ous that the symmetry between the ciu-ve A = 1.0 and 
A = —1.0 is due to the symmetry of geometry. Now we 
will focus on the case of A = —1.0, namely, the right 
side is steeper than the left one. In the adiabatic limit 
w — >■ 0, the vibration can be expressed by two opposite 
static driving — &o and bo and the nonequilibrium sys- 
tem reduces to the system with two equilibrium states. 



namely, v = 6o) + ''^(^o)] = 0. For low frequencies, 

the particles get enough time to reach the whole area and 
the wall will act on the particles. In this case, the proba- 
bility of the particles acting on the slanted wall (the left 
wall) is larger than that on the steeper wall (the right 
wall), resulting in a positive current. On increasing the 
frequency w. due to high frequency, the Brownian parti- 
cles do not get enough time to reach the slanted wall (The 
area of the slanted side is larger than that of the steep 
one) . Therefore, the probability of the particles acting on 
the slanted wall is smaller than that on the steep wall, 
yielding a negative current. When the vibration changes 
very fast w — >■ oo, the particles can not experience the 
changes of the wall and the system reduces to a equilib- 
rium system, so the current goes to zero. Therefore, the 
direction of the current can be controlled by changing the 
vibrating frequencies. 

In order to give more detail insight into the character- 
istic frequencies ujl, ujc, and LUfj, we corroborate the nu- 
merical estimates of the frequencies for a larger range of 
parameter values. The numerical estimates are shown in 
Fig 4. We find that the characteristic frequencies increase 
with temperature T and the length of the period L and 
decrease with the vibrating amplitude bo . The estimated 
exponential relations are also shown in the figure. Prom 
the figure we can obtain the approximate scaling results 
for UJL oc rO-646-i-57i^-o.95^ T°-'%'-^' L-^'^^ , and 

u)R oc T°-^25-i.072^-i.4i jg obvious that these char- 
acteristic frequencies will shift to large values for high 
temperature and to small values for large amplitude of 
the vibration. 

In order to investigate how the temperature affects the 
transport, we also show the mean velocity as a function 
of temperature in Fig. 5(a) at A = —1.0 and ui = 10.0. 
The diffusion from the temperature has two roles: jump- 
ing from one cell to another in longitudinal direction and 
making the particles acting on the wall in transverse di- 
rection. When T — ^ 0, the particles cannot reach the wall 
and the effect of wall disappears and there is no current. 
When r — > oo, the effect of the wall vibration disappears 
and the current goes to zero, also. Therefore, there is 
an optimized value of T at which the mean velocity v 
takes its maximum value, which indicates that the ther- 
mal noise may facilitate the particles transport. In order 
to give a feeling of the generality of the optimized value of 
temperature, we also study the dependence of Tc on the 
vibrating frequency and amplitude. Prom Fig. 5 (b) and 
(c), we can obtain the approximate scaling results for the 
optimized temperature Tc oc boOJ°'^^. The value of the 
optimized temperature will become larger for high driv- 
ing frequency. Remarkably, the optimized temperature 
is independent of the driving amplitude. 

Figure 6 (a) shows the mean velocity t; as a function 
of the controlling parameter b. The parameter b deter- 
mines the radius at the bottleneck. The tube may be 
close for too small values of b. If the radius at the bot- 
tleneck is small few particles can pass from one cell to 
another, so the mean velocity is small. When the bottle- 
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FIG. 4: Dependence of characteristic frequencies ujl, (jJc, and 
LOR on the different parameters. (a)On temperature T at a = 
0.5, b = 0.9, bo = 0.3, A = -1.0, and L = 1.0. (b)On the 
vibrating ampHtude 6o at T = 0.05, a = 0.5, b = 0.9, A = 
-1.0, and L = 1.0. (c)On the period L at T = 0.05, a = 0.5, 
b = 0.9,6o = 0.3, and A = —1.0. The solid Unes are the fitting 
Unes. 



FIG. 5: (a)The mean velocity tj as a function of temperature 
T at a = 0.5, b = 0.9, bo = 0.3, A = -1.0, cj = 10.0, and 
L = 1.0. Tc is the optimized value of temperature for the 
peak. (b)Dependence of Tc on the driving frequency w at 
a = 0.5, b = 0.9, bo = 0.3, and A = —1.0. (c)Dependence of 
Tc on the driving amplitude 6o at a = 0.5, b = 0.9, w = 10.0, 
and A = —1.0. The solid Unes are the fitting lines. 
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(b) 




neck has infinite radius (b — > oo), the tube reduces to a 
straight one and the effect of tube shape disappears, so 
the current goes to zero. Therefore, there exists a value 
of be at which the mean velocity takes its maximum. 
We also study how the temperature T and the vibrating 
frequency u affect the optimized value be- Prom the nu- 
merical simulations, we also find the approximate scaling 
results for be oc T~^-^^oj^-^^ . It is found that be decreases 
with temperature T and increases wit the vibrating fre- 
quency ui. 
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FIG. 6: (a)Thc mean velocity t; as a function of parameter 
b at a = 0.5, bo = 0.3, A = -1.0, uj = 10.0, T = 0.05, and 
L — 1.0. be denotes the optimized value of b for the peak. 
(b)Dependence of be on temperature T at a = 0.5, 6o = 0.3, 
A = -1.0, cj = 10.0, and L = 1.0. (b)Dependence of be 
on the driving frequency u) at a = 0.5, bo = 0.3, A = —1.0, 
T = 0.05, and L = 1.0. The solid lines are the fitting lines. 



FIG. 7: The mean velocity « as a function of vibrating am- 
plitude 6o at o = 0.5,6 = 0.9, A = -1.0, w = 10.0, T = 0.05, 
and L = 1.0. 

Figure 7 gives the mean velocity w as a function of the 
vibrating amplitude 6o- When 6o — >^ 0, the particles can 
not experience the wall vibrating, so the current goes 
to zero. On increasing the vibrating amplitude b{), the 
mean velocity increases. However, for large values of bo, 
the tube may be close at some time, then the current 
decreases. We must point out that for very large value 
of 6o) the profiles of tube wall may be exchanged, this is 
not true for a real system. 

IV. CONCLUDING REMARKS 

In this paper, we study the transport of overdamped 
Brownian particles moving in an asymmetrically periodic 
tube with the wall vibration. From the ntimerical sim- 
ulation we can find that the transverse wall vibration 
can induce a net longitudinal current in an asymmetri- 
cally periodic tube. We can have current reversals by 
changing the sign of A, the asymmetry of the tube. The 
sign of mean velocity for low frequencies is opposite to 
that for high frequencies. Therefore, we can change the 
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direction of the current by suitably tailoring the vibrat- 
ing frequency. Thus, varying the asymmetry of the tube 
and the vibrating frequency are the two ways of inducing 
current reversals. These peculiar reversals are due to dif- 
ferent strength and asymmetry of relaxation processes. 
The symmetry of relaxation processes changes with the 
system parameters, so the current may changes its di- 
rection when the parameters are changed. In addition, 
there exists an optimized value of temperature at which 
the mean velocity takes its maximum value, which indi- 
cates that the thermal noise may facilitate the particles 
transport. 

Clearly, the model is too simple to provide a real- 
istic description of real systems, however, the results 
we have presented have a wide application in may pro- 



cesses, such as diffusion of ions and macromolecular so- 
lutes through the channels in biological membranes^, 
transport in zeolites^^and nanostructures of complex 
geometry^-, controlled drug releasees, and diffusion in 
man-made periodic porous materiala^^. It is very impor- 
tant to understand the novel properties of these confined 
geometries, zeolites, biological channels, nanoporous ma- 
terials, and microfluidic devices, as well as the transport 
behavior of species in these systems. 
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